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Abstract 

Arnold showed that the Euler equations of an ideal fluid describe 
geodesies on the Lie algebra of incompressible vector fields. We gener- 
alize this to fluids with dissipation and Gaussian random forcing. The 
dynamics is determined by the structure constants of a Lie algebra, along 
with inner products defining kinetic energy, Ohmic dissipation and the 
covariance of the forces. This allows us to construct tractable toy models 
for fluid mechanics with a finite number of degrees of freedom. We solve 
one of them to show how symmetries can be broken spontaneously.In 
another direction, we derive a deterministic equation that describes the 
most likely path connecting two points in the phase space of a randomly 
forced system: this is a WKB approximation to the Fokker-Plank-Kramer 
equation, analogous to the instantons of quantum theory. Applied to hy- 
drodynamics, we derive a PDE system for Navier-Stokes instantons. 



1 The Euler Equations 

The velocity of an ideal fluid satisfies 

dm k 



= -v k 3 k v t - V iP (1) 
along with the auxiliary condition of incompressibility 

d k v k = 0. (2) 

Arnold 1 advocated an elegant geometric interpretation for these equations: 
they describe geodesies on the Lie algebra of incompressible vector fields. 
Recall that the commutator 



[u,v 



k = u 3 d 3 v k -v 3 d jU k . (3) 



of two incompressible (i.e., zero divergence) vector fields is still incompressible: 
such vector fields form a Lie algebra. Now, any inner product on a Lie algebra 
defines a left-invariant Riemannian metric on the corresponding Lie group [2]. 
If the Lie algebra is non-abelian, this is a curved metric. But, the evolution of 
the tangent vectors (momenta) along a geodesic decouple from the co-ordinate 
variables on the group, because of the left invariance under the group action. 
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1.1 Hamiltonian Formalism 



A more physical point of view is to think of the Lie algebras as defining Poisson 
brackets on the (angular) velocities 

{v a ,V b } = fa b Vc (4) 

The inner product defines the kinetic energy 

E = ^G ab v a v b (5) 

The resulting Hamilton's equations are the geodesic equations of the underlying 
left-invariant metric on the Lie group: 

^ = { Va ,E} = ti b G bd v c v d (6) 

More precisely, this is how the tangent vector evolves along a geodesic. 
For an ideal fluid, the kinetic energy provides a natural inner product: 

1 r. 



E=-J v'dx (7) 

Arnold's observation is that the Euler equations (fTJ) are precisely the geodesic 
equations defined by this metric ([7]) on the Lie algebra above. A more 
detailed derivation can be found in (3J |4] - 

More generally, the kinetic energy is a quadratic function of the velocities, 
so defines a Riemannian metric on the configuration space. In the absence 
of potential energy, the system moves along the geodesies of this metric: the 
time variable is proportional to the arc-length. Arnold's point is that an ideal 
fluid is a conservative (hamiltonian) system without any potential energy; so 
it moves along the geodesies of the metric defined by its kinetic energy. A 
compressible fluid has potential energy, stored in the deformation of density 
from its equilibrium value; so its time evolution is not along geodesies. Also, a 
viscous fluid is not conservative, so the Navier-Stokes equations are not geodesic 
equations either. 

The curvature of a Riemannian manifold describes the stability of geodesic 
flows [SJ If the sectional curvature R(u,v,u,v) < two geodesies that start 
at the same point along tangent vectors u, v will diverge from each other. Arnold 
showed that in all except a finite number of directions, the sectional curvature 
is negative for the Euler flow, perhaps explaining the instability of fluid flow. 



1.2 Example: Euler Equations of A Rigid Body 

The most well known system to which we can apply the hamiltonian formalism 
above is the rigid body. The Lie algebra is that of rotations 

{Li, Lj} = eijkLk (8) 

where Li are the components of angular momentum as measured in the comoving 
frame of the body. Because this is not an inertial frame, the angular momentum 
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may not be conserved even in the absence of external torque. The kinetic energy 
is (in the frame where the moment of inertia tensor is diagonal) 

l - [GtLl + G 2 L 2 2 + G 3 Lf\ (9) 

leading to the Euler equations of a rigid body. These are equations for geodesies 
on a triaxial ellipsoid: the group 5*0(3) with a left-invariant metric. 
The curvatures in the principal planes are 

K23 = 4G^G 3 (10) 

and cyclic permutations thereof. 



1.2.1 An instability of coin flip 

For a cylinder with 




(e.g., coin) rotations with axis in the plane of symmetry are unstable: curva- 
ture is negative. To see this instability we would have to mark a point of the 
circumference of a coin. The rotation of a coin around an axis that lies on the 
plane of the coin is unstable (i.e., the marked point will deviate) with respect to 
a perturbation of the axis to another one which also lies in the plane of the coin. 
An experiment to see this phenomenon is not too hard to set up with modern 
equipment (video cameras) but appears not to have been performed ever. 



1.3 An Anomaly 

Recall that in hamiltonian mechanics there is always a volume form in phase 
space that is invariant under time evolution: this is Liouville's theorem. In our 
case we can ignore the position variables as they decouple from the evolution of 
velocities. If we simply take the divergence, 

- fa (12) 

dv a ~ U dv b (12) 

If the Lie algebra is unimodular, 

r ab = o, (13) 

^ = 0. (14) 
ov a 

Thus the constant measure is invariant. Any semi-simple Lie algebra is unimod- 
ular. But some nilpotent and solvable ones are not. 

Finite dimensional non-unimodular algebras have some other measure /i on 
momentum space that is preserved by the geodesic flow: 

« = 0. (15, 



this is zero 
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We will find it explicitly in some special cases below. This subtle point, that the 
constant measure in velocity space may not be invariant under time evolution, is 
reminiscent of an 'anomaly' in quantum field theory. The anomalies of quantum 
field theory are determined by the cohomology of the Lie algebra twisted by 
its action on some space of fields; the invariant measure \x also determines an 
element of the (untwisted9 Lie algebra cohomology; when this element is not 
non-trivial we have an anomaly as well: /i is not constant. We will see that, as 
in field theory, such an anomaly can lead to spontaneous symmetry breaking. 

1.4 Example: The Poincare Half Plane 

The only non-abelian Lie algebra in two dimensions is the affine Lie Algebra 

{uo,Ui} = Ui- (16) 

This algebra is not unimodular. The invariant measure on velocity space is 

a dv dv 1 
ixav = . (17) 

Vl 

The point is that vq generates scale transformations in vi , so to be invariant we 
must divide the naive measure dv^dvi by v\. This example therefore is a toy 
model where we can investigate the consequences of the anomaly we mentioned 
earlier. 

Any quadratic form on this Lie algebra can be reduced to 

E=\[vl+vi]. (18) 

The corresponding group manifold is the upper half-plane, with the Poincarcmctric. 
The geodesic equations are, using the above hamiltonian formalism 

-df = - v ^^r = VoVl (19) 

Solutions are semi-circles in momentum space: 

t>o = — ptanhpi, v\ — psech pt. (20) 
We will see later how dissipation and random forcing affect this dynamics. 



2 Geometry of Dissipative Motion 

A conservative dynamical system is described by a Poisson bracket and hamil- 
tonian (energy). Dissipation is described by (a possibly degenerate, but always 
non-negative) tensor T ab , which determines the gradient of energy: 

l-««- r *i <21) 

It is straightforward to check that the Energy decreases: 

dE dEdE 
dt d^ a <9£ fc 
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In the cases we study in this paper, the Poisson bracket is linear ( a Lie algebra), 
the energy is a quadratic function and r a f, is constant. 

^ = -r afc G bc i; c + f ab G bd v c v d (23) 

Even this special case is quite complex: it includes Navier- Stokes, with the 
choices 

E=\{v,Gv) = ^J v 2 dx, {v, Tv) = 1 J {diVjf dx. (24) 

Thus the quadratic form of kinetic energy is the L 2 norm while that of 
dissipation is the H 1 Sobolev norm. 



3 The Langevin And Fokker-Planck-Kramers 
Equations 

A standard model of randomness [B] is a Gaussian white noise: a force whose 
correlations vanish fast compared to the dynamical and dissipative time scales. 
If the random force is a sum of a large number of more or less independent 
forces, it will be a Gaussian. The covariance of the Gaussian introduces a third 
positive tensor into the game. 

^ = -T ab G bc v c + f ab G bd v c v d + Va , < Va(t)vb(t') >= D ab 5(t - t'). (25) 

This stochastic differential equation (Langevin equation) implies the Fokker- 
Planck-Kramers Equation for the evolution of the probability density of velocity: 



dP _ d 
^ dt dv a 



H |^ab|£ - {-T ab G bc v c + f r ab G bd v c v d ) P I 



(26) 



Here, /i is the invariant density of the ideal dynamics. 

The traditional derivation [7 4 [5] of this FPK equation does not include the 
density ji : it is implicitly assumed that the invariant (Liouville) density is the 
standard one. In order that the total probability be conserved: 

- / fiPdv = 0. (27) 



dt . 

we must order the differentiations and insert factors of /i as above. This is one 
of our innovations in this paper. 

3.1 Maxwell-Boltzmann is Equilibrium Solution for Uni- 
modular Lie Algebras 

If the Einstein relation 

PD ab = T ab (28) 

holds, P = e~P E is a static solution of the FPK equations. But it may not be 
normalizable, because J \ie~@ E dv might diverge. For unimodular Lie algebras, 
H = 1 and the Maxwell-Boltzmann distribution P = e~3^( 1, > Gt ') is an equilibrium 
solution. An example is the random motion of a rigid body. This case has been 
studied in detail by astronomers interested in the motion of dust grains in the 
interstellar medium [5]. We now turn to an example with /i =/= 1. 
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3.2 The Langevin Equation on the Half-Plane 



The geodesic equations still follow for zero dissipation if we choose as hamilto- 
nian any function H(p) of p = yv^ +v\. We can choose H(p) such that the 
FKP equation is easier to solve. But physically, the kinetic energy would still 
The Fokkcr-Planck-Kramcr reduces to 



be E = \p 2 



dQ 

at 



D 



8 2 Q 
dp 2 



dQ 

19 dp + - 2 *° 



D d 



cosh 6— 
dd 



Q (29) 



This equation is separable for H — — -. In fact the solution can be found in 
closed form in terms of spheroidal functions. 

Numerical solution of toy model shows that the details of the dynamics 
are not important: the vestige of the dynamics that is important to the long 
term behavior is the invariant measure /i. We can even set H = and get a 
qualitatively similar equilibrium distribution. In this case the solution is quite 
simple [TU]: 



P(vo,vi) 



dvndvi Vi _j3, 

= -~e 2' 

vi p z 



1 - J-e^erf ( ^l) 
^p G er 'l 2 j 



dv^dvi 



(30) 



The main point is that this solution does not have a peak at minimum energy 
vq = vi = 0; instead the peak is at some finite value of vi. The anomaly drives 
a spontaneous breaking of the rotation invariance in the energy E — \ [vq + vf]. 
Even if the two components of velocity appear symmetrically in energy, their 
commutation relations are different. These commutation relation determine the 
measure of integration which is singular at vi = 0. It is this singularity that 
pushes the most likely equilibrium value of Vi away from zero. 

Thus we see that anomaly we discovered can drive a spontaneous symmetry 
breakdown: a steady state flow along the xi axis is generated in the equilib- 
rium. Perhaps this is a toy model for the formation of flows such as hurricanes 
and ocean currents which persist for long times despite the instability of the 
underlying fluid system. 



4 The WKB Approximation to the FPK Equa- 
tion 

It is often important to know the probability that a certain initial condition 
will evolve to a given final condition under the influence of random forces. It 
might be important to compute this even when this probability is small: there 
could be final states that are catastrophic so that even a small chance of them 
happening cannot be ignored. If the probability is small, we can use the WKB 
approximation as for the Schrodinger equation; put P = e~* and ignore second 
derivatives of $: 

a + D -^^l + H 1 - ** + ^ G = < 31) 

These are the Hamilton- Jacobi equations of the hamiltonian 

H = D ab w a w b + (-T ab G bc v c + f c ab G bd v c v d ) w a , (32) 
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where w a are the canonical conjugates of the variables v a '- 

{v a ,v b } 1 = = {w a ,w b } , {v a w b } 1 = S b (33) 

Kramer used this WKB approximation to the FPK in chemical kinetics; we 
extend them to a more general framework. In quantum mechanics, the WKB 
approximation works when the quantum tunneling amplitudes are small; the 
path with the greatest probability is called the instanton. The analogue in our 
case is the solution to the Hamilton's equations following from the above efectie 
hamiltonian. 

(11 ) (11 11^ 

= 2D ab w b -T ab G hc v c + f c ah G bd v c v d , — = (T ab G bc ~ 2f c ab G bd v d ) w a (34) 

Notice that this effective hamiltonian mechanics of rare transitions has twice 
the number of degrees of freedom as the original conservative dynamics we 
started with: we can specify initial and final conditions independently. The 
action of the path that solves this equation with given initial and final states 
gives the log of the transition probability: 

$(v',v,T)= [ [w a v a - H(v,w)]dt (35) 
Jo 

4.1 Instantons of Random Motion on the Half-Plane 

With an isotropic dissipation tensor, these instanton equations are quite simple 
for the half-plane: 

dv 2 dvi 

— =w ~ 7^0 - v x , — =wi- 7«i + v vi (36) 
dw dw-i 

—— = j Wq - viwi, —— = 71W1 + 2viw - v wi (37) 
dt dt 

We could not find the general solution to this system. But notice the special 

solution with w a = 2"/v a . This is just the time reverse of the dissipative solution! 

This solution leads to the Maxwell-Boltzmann distribution: agrees with the 

exact solution for large energies. WKB is not good for small energies, where 

the exact answer departs from MB distribution. This is the phenomenon of 

spontaneous symmetry breakdown we discovered in the last section. 



4.2 Navier-Stokes Instantons 

We can now derive the PDE describing the most likely time evolution of a 
randomly forced dissipative fluid: 

=Fi + 7<9 2 i; l - VjdjVi - dip, d k v k = 0. (38) 

—1 -jd 2 Wi + wjdjVi + vjdjWi di(f), d k w k = 0. (39) 
The covariance of the fluctuations is likely to be an integral operator [6] 



Jl-Kik-(x — x 1 ) 



kikj 

~fc2~ 



D(k)dk. 



(40) 
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If we choose 



D(k) = \k\ 



-2 



(41) 



we get an auxiliary differential equation to complete the above system 



With this choice, the random forcing grows with spatial scale, while dissipation 
is more important at short scales: what is expected in fluid mechanics. The 
above equation should be a good approximation to calculate probabilities and 
paths of large deviations from Navier-Stokes. Such rare events are important 
in many situations: the path of a hurricane, or the chances of a tsunami or 
a rogue wave. Similar instanton equations can be derived for Rotating Shal- 
low Water equations, the Stochastic Loewner Evolution [TT], two dimensional 
incompressible flow etc. We hope to return to these questions in later publica- 
tions. Numerical solution of our system is not much harder than the solution of 
Navier-Stokes itself: only twice as many variables; moreover, the new variables 
appear linearly. 
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d 2 F l = 2Wi. 



(42) 
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